grad(sin(pi*q(x,y,z))*sinh(pi*q(x,y,z))^(-3/2))